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Transverse spin structure of the nucleon from lattice QCD simulations 
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We present the first calculation in lattice QCD of the lowest two moments of transverse spin 
densities of quarks in the nucleon. They encode correlations between quark spin and orbital angular 
momentum. Our dynamical simulations are based on two flavors of clover-improved Wilson fermions 
and Wilson gluons. We find significant contributions from certain quark helicity flip generalized 
parton distributions, leading to strongly distorted densities of transversely polarized quarks in the 
nucleon. In particular, based on our results and recent arguments by Burkardt [Phys. Rev. D 72 
(2005) 094020] , we predict that the Boer- Mulders- function hi , describing correlations of transverse 
quark spin and intrinsic transverse momentum of quarks, is large and negative for both up and down 
quarks. 



Introduction. — The transverse spin (transversity) 
structure of the nucleon received a lot of attention in 
recent years from both theory and experiment as it pro- 
vides a new perspective on hadron structure and QCD 
evolution (for a review see [l|). A central object of inter- 
est is the quark transversity distribution Sq(x) — hi(x), 
which describes the probability to find a transversely po- 
larized quark with longitudinal momentum fraction x in 
a transversely polarized nucleon [2]. Much progress has 
been made in the understanding of so-called transverse 
momentum dependent PDFs (tmdPDFs) like e.g. the 
Sivers function /f^x, fcj_) |3|, which measures the cor- 
relation of the intrinsic quark transverse momentum k± 
and the transverse nucleon spin S±, as well as the Boer- 
Mulders function h^(x, kj_) [J], describing the correlation 
of k±_ and the transverse quark spin sj_. While the Sivers 
function begins to be understood, still very little is known 
about the sign and size of the Boer-Mulders function. 



experimentally accessible asymmetries. 

Lattice calculations give access to re-moments of trans- 
verse quark spin densities [B[ 
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where m is the nucleon mass. The ^-dependent nucleon 
generalized form factors (GFFs) A n o(b\), AT n o{b\), . . . 
in Eq. (Q} are related to GFFs in momentum space 
A n o(t), Axnoit), ... by a Fourier transformation 



A particularly promising approach is based on 
3-dimensional densities of quarks in the nucleon, 
p{x, b±, sj_, S±) |5j, representing the probability to find a 
quark with momentum fraction x and transverse spin s± 
at distance b± from the center-of-momentum of the nu- 
cleon with transverse spin S± . As we will see below, these 
transverse spin densities show intriguing correlations of 
transverse coordinate and spin degrees of freedom. Ac- 
cording to Burkardt @, 0], they are directly related to 
the above mentioned Sivers- and Boer-Mulders- functions. 
Our lattice results on transverse spin densities therefore 
provide for the first time quantitative predictions for the 
signs and sizes of these tmdPDFs and the corresponding 
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where Aj^ is the transverse momentum transfer to the 
nucleon. Their derivatives are defined by /' = d& f and 
Ab ± f = 49 fc ^ (b\d b 2)f. The generalized form factors in 
this work are directly related to x-moments of the cor- 
responding vector and tensor generalized parton distri- 
butions (GPDs) (for a review see 3). The probability 
interpretation of GPDs in impact parameter space has 
been first noted in [9( . Apart from the orbitally symmet- 
ric monopole terms in the second line of Eq.([T]), there are 
two dipole structures present in the third line of Eq.(p}, 
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FIG. 1: Results for the generalized form factors -B T ( n = l!2 )o(£). 
The corresponding p-pole parametrizations are shown by the 
shaded error bands. 



31 s\ and 6j_e Jl 5j_. The fourth line in Eq.fT]) corre- 
sponds to a quadrupole term. The (derivatives of the) 
three GFFs B n0 (b±), BT n o(b±) and AxnoipA.) thus de- 
termine how strongly the orbital symmetry in the trans- 
verse plane is distorted by the dipole and the quadrupole 
terms. 

The GFFs A n o(t), Axnoit), . . . parametrize off- forward 
nucleon matrix elements of certain local quark operators. 
For the lowest moment n = 1 one finds Aio(t) = F±(t), 
B 10 (t) = F 2 (t) and A TW (t) = g T (t) where Fx, F 2 and 
gx are the Dirac, Pauli and tensor nucleon form factors, 
respectively. A concrete example of the corresponding 
parametrization for n = 1 is given by 1^, 11 1 
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where O^f = 



qa^ v ^ 5 q is the lowest element of the tower 
of local leading twist tensor (quark helicity flip) oper- 
ators. Parametrizations for higher moments n > 1 in 
terms of tensor GFFs and their relation to GPDs are 
given in ll|. As it is very challenging to access tensor 
GPDs in experiment [12], input from lattice QCD calcu- 
lations is crucial in this case. 

Simulation results. — Our lattice calculations are based 
on configurations generated with nj — 2 dynamical non- 
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FIG. 2: Study of discretization errors of the tensor charge 
j4tio(£=0) = gT(t=Q) for up- and down-quarks at a pion mass 
of m n w 600 MeV. 



perturbatively 0{a) improved Wilson fermions and Wil- 
son gluons. Simulations have been performed at four 
different couplings /3 = 5.20, 5.25, 5.29, 5.40 with up 
to five different k = K sca values per (3, on lattices of 
V x T = 16 3 x 32 and 24 3 x 48. The lattice spacings 
are below 0.1 fm, the range of pion masses extends down 
to 400 MeV and the spatial volumes are as large as 
(2.1 fm) 3 . The lattice scale a in physical units has been 
set using a Sommer scale of ro — 0.467 fm 1J, The 
computationally demanding disconnected contributions 
are not included. We ex pec t, however, that they are 
small for the tensor GFFs [15j . We use non-perturbative 
renormalization [l6| to transform the lattice results to 
the MS scheme at a scale of 4 GcV 2 . The calculation 
of GFFs in lattice QCD follows standard methods (see, 

e.g., EHEi). 

In Fig. [TJ we show as an example results for the GFFs 

^T(n=i,2)oW' corresponding to the lowest two moments 

n = 1,2 of the GPD ~E^ d (x,i,t) [2(J, as a function of 
the momentum transfer squared t, for a pion mass of 
m^ rs 600 MeV, a lattice spacing of a m 0.08 fm and 
a volume of V w (2 fm) 3 . For the extrapolation to the 
forward limit (t — 0) and in order to get a functional 
parametrization of the lattice results, we fit all GFFs us- 
ing a p-pole ansatz F(t) = Fq/(1 — (t/(pm^)) p with the 
three parameters Fq = F(t=0), m p and p for each GFF. 
We consider this ansatz [2l| to be more physical than 

oc m" 1 is inde- 



previous ones as the rms-radius (r 2 ) 1 / 2 
pendent of p. It turns out that in most cases the statistics 
is not sufficient to determine all three parameters from 
a single fit to the lattice data. For a given generalized 
form factor, we therefore fix the power p first, guided by 
fits to selected datasets, and subsequently determine the 
forward value Fq and the p-pole mass m p by a full fit to 
the lattice data. Some GFFs show a quark flavor depen- 
dence of the value of p, which has already been observed 
in [22| for the Dirac form factor. For the examples in 
Fig. [U we find for u-quarks B^ w (t=0) = 3.34(8) with 
to p = 0.907(75) GeV, Bl 20 {t=0) = 0.750(32) with m p = 

1.261(40) GeV and for d-quarks B^ w (t=0) = 2.06(6) 
with m p = 0.889(48) GeV, B^ 20 (t=0) = 0.473(22) with 
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FIG. 3: Pion mass dependence of the generalized form factors 
B T ( n=1 ,2)0 (* = 0) f° r up-quarks. The shaded error bands show 
extrapolations to the physical pion mass based on an ansatz 
linear in m%. The symbols are as in Fig. 2. 



to p = 1.233(27) GeV (all for p = 2.5). We have checked 
that the final p-pole parametrizations only show a mild 
dependence on the value of p chosen prior to the fit. In 
order to see to what extent our calculation is affected 
by discretization errors, we plot as an example in Fig. [2] 
the tensor charge ATio{t=0) = gT(t=0) versus the lat- 
tice spacing squared, for a fixed w 600 MeV. The 
discretization errors seem to be smaller than the statis- 
tical errors, and we will neglect any dependence of the 
GFFs on a in the following. Taking our investigations of 
the volume dependence of the nucleon mass and the axial 
vector form factor ID, 2^ as a guide, we estimate that 
the finite volume effects for the lattices and observables 
studied in this work are small and may be neglected. 

As an example of the pion mass dependence of our 
results, we show in Fig. [3]the GFFs Bt(h=i,2)o(' = ") ver " 
sus m 2 . Unfortunately we cannot expect chiral pertur- 



bation theory predictions [2J| to be applicable to most 
of our lattice data points, for which the pion mass is 
still rather large. To get an estimate of the GFFs 
at the physical point, we extrapolate the forward mo- 
ments and the p-pole masses using an ansatz linear in 
77i 2 . The results of the corresponding fits are shown as 
shaded error bands in Fig. El At mf ys = 140 MeV, 



we find B T10 (t=0) 



B T2Q (t-0) 



2.93(13) 
0.420(31), B d T2Q (t=0) 



B^w (*=0) 



1.90(9) and 
0.260(23). These 



comparatively large values already indicate a significant 
impact of this tensor GFF on the transverse spin struc- 
ture of the nucleon, as will be discussed below. Since the 
(tensor) GPD Et can be seen as the analogue of the (vec- 
tor) GPD E, we may define an anomalous tensor mag- 
netic moment • k / J dxET{x,^,t=Q) = Byio(t=0), 
similar to the standard anomalous magnetic moment 
k = J dxE(x,£,t=0) = B w (t=0) = F 2 (t=0). While the 
u- and d-quark contributions to the anomalous magnetic 
moment are both large and of opposite sign, 
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FIG. 4: Lowest moment (n = 1) of the densities of un- 
polarized quarks in a transversely polarized nucleon (left) 
and transversely polarized quarks in an unpolarized nucleon 
(right) for up (upper plots) and down (lower plots) quarks. 
The quark spins (inner arrows) and nucleon spins (outer ar- 
rows) are oriented in the transverse plane as indicated. 
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the anomalous tensor magnetic moment for both flavors, 



values have been obtained in a recent model calculation 
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[25J|. Large N c considerations predict k t p ~ n, T 

Let us now discuss our results for p n {b±, s±, S±) in 
Eq. ([1]). For the numerical evaluation we Fourier trans- 
form the p-pole parametrization to impact parameter 
space. The parametrizations of the impact param- 
eter dependent GFFs then depend only on the p-pole 
masses m p and the forward values Fq. Before showing 
our final results, we would like to note that the mo- 
ments of the transverse spin density can be written as 
sum/difference of the corresponding moments for quarks 
and antiquarks, p n = p™ + (— l)™p^, because vector and 
tensor operators transform identically under charge con- 
jugation. Although we expect contributions from anti- 
quarks to be small in general, only the n-even moments 
must be strictly positive. In Fig. 01 we show the lowest 
moment n = 1 of spin densities for up and down quarks 
in the nucleon. Due to the large anomalous magnetic 
moments k u ' , we find strong distortions for unpolarized 
quarks in transversely polarized nucleons (lefVpart of the 
figure). This has already been discussed in [6j, and can 
serve as a dynamical explanation of the experimentally 
observed Sivers-effect. Remarkably, we find even stronger 
distortions for transversely polarized quarks s± — (s x , 0) 
in an unpolarized nucleon, as can be seen on the right 
hand side of Fig. [4j The densities for up and for down 
quarks in this case are both deformed in positive b y direc- 
tion due to the large positive values for the tensor GFFs 
B^ 1Q (t=0) and B rio (f=0), in strong contrast to the dis- 
tortions one finds for unpolarized quarks in a transversely 
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FIG. 5: Second moment (n = 2) of transverse spin densities. 
For details see caption of Fig. 4. 



Conclusions. — We have presented first lattice results 
for the lowest two moments of transverse spin densities 
of quarks in the nucleon. Due to the large and positive 
contributions from the tensor GFF Bq^no for up and for 
down quarks, we find strongly distorted spin densities for 
transversely polarized quarks in an unpolarized nucleon. 
According to Burkardt [7J , this leads to theprediction of 
a sizable negative Boer-Mulders function [4j for up and 
down quarks, which may be confirmed in experiments at 
e.g. JLab and GSI/FAIR 

The numerical calculations have been performed on 
the Hitachi SR8000 at LRZ (Munich), the apeNEXT 
at NIC/DESY (Zeuthen) and the BlueGene/L at 
NIC/FZJ (Jiilich), EPCC (Edinburgh) and KEK (by 
the Kanazawa group as part of the DIK research pro- 
gramme). This work was supported by DFG (Forscher- 
gruppe Gitter-Hadronen-Phanomenologie and Emmy- 
Noether programme), HGF (contract No. VH-NG-004) 
and EU I3HP (contract No. RII3-CT-2004-506078). 



polarized nucleon. All of these observations are actually 
quite plausible, because there is no gluon transversity 
which could mix with quarks under evolution. Therefore, 
the transverse spin structure is much more valence-like 
than the longitudinal one, which is strongly affected by 
the negative sign of the photon-gluon contribution. Thus 
the transverse quark spin and the transverse quark or- 
bital angular momentum simply seem to be aligned. The 
fact that the u-quark (d-quark) spin is predominantly 
oriented parallel (antiparallel) to the nucleon spin then 
explains why densities of quarks with spin in ir-direction 
in an unpolarized nucleon (moving towards the observer 
in z-direction) are larger in the upper half plane. How- 
ever, the contributions from down quarks with spin in 
(— x)-direction dominate in a polarized nucleon with spin 
along the (+a;)-axis, such that the orbital motion around 
the (— x)-direction leads to a larger down quark density 
in the lower half plane. It has been argued by Burkardt 
that the deformed densities on the right hand side 
of Fig. [4] are related to a non- vanishing Boer-Mulders 
function [4| hj~ which describes the correlation of in- 
trinsic quark transverse momentum and the transverse 
quark spin sj_. According to 0] we have in particular 
KT —hi. If this conjecture is correct our results imply 
that the Boer-Mulders function is large and negative both 
for up and down quarks. The fact that the correlation of 
quark and nucleon spin is not 100 percent explains why 
the deformation is more pronounced in the Boer-Mulders 
than in the Sivers case. 

Fig.[5]shows the n = 2-moment of the densities. Obvi- 
ously, the pattern is very similar to that in Fig. [4j which 
supports our simple interpretation. The main difference 
is that the densities for the higher n — 2-moment are 
more peaked around the origin = as already ob- 
served in [271 for the vector and axial vector GFFs. 
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